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Abstract To date, infectious bronchitis vims (IBV) is po¬ 
tentially found in wild birds of different species. This 
work reports the survey of coronavimses in wild birds 
from Madagascar based on the targeting of a conserved 
genome sequence among different groups of CoVs. Phy¬ 
logenetic analyses revealed the presence of 
gammacoronaviruses in different species of Gmiformes, 
Passeriformes, Ciconiiformes, Anseriformes, and 
Charadriiformes. Furthermore, some sequences were relat¬ 
ed to various IBV strains. Aquatic and migratory birds 
may play an important role in the maintenance and spread 
of coronavimses in nature, highlighting their possible con¬ 
tribution in the emergence of new coronavims diseases in 
wild and domestic birds. 
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Introduction 

Coronavimses (CoVs) are enveloped vimses within the 
Coronaviridae family with a positive-sense RNA genome 
ranging from 27 to 30 kb (Cavanagh 1997). They are divided 
into four genera: Alpha-, Beta-, Gamma- and 
Deltacoronavirus. 

Alpha- and betacoronavimses have been isolated from 
mammals; gammacoronaviruses , being the most representa¬ 
tive member infectious bronchitis vims (IBV), have also been 
detected in the beluga whale (Mihindukulasuriya et al. 2008), 
the Asian leopard cat (Dong et al. 2007), and in several species 
of wild birds (Chu et al. 2011; Hughes et al. 2009; Jonassen 
et al. 2005; Muradrasoli et al. 2010; Woo et al. 2009). Vimses 
of the new proposed genus Deltacoronavirus were found pri¬ 
marily in mammals, terrestrial and aquatic birds (Chu et al. 
2011; Dong et al. 2007; Woo et al. 2009). 

In animals, Co Vs are associated to respiratory and intesti¬ 
nal infections, as well as hepatic, renal, and neurological dis¬ 
orders (Cavanagh 2005). Human coronavimses infections 
were generally related to the common cold, but the emergence 
and pandemic potential of severe acute respiratory syndrome 
(SARS) in 2003 and the Middle East respiratory syndrome 
(MERS) in 2012, both belonging to genus Betacoronavirus , 
led to the discovery of several Co Vs hosted by humans and 
animals (Hilgenfeld and Peiris 2013). 

Over the past years, the wildlife has been under epidemio¬ 
logical surveillance worldwide, as they can play an important 
role as reservoir of emerging vimses that may pose a risk to 
mankind and threaten wildlife itself. Wild birds are known to 
be important reservoirs of avian influenza A vims (Brown and 
Stallknecht 2008; Spackman 2009), as well as other respira¬ 
tory and enteric vimses like IBV (Cavanagh and Gelb 2008). 
In poultry, the disease is controlled through the use of vaccines 
and biosafety measures, including avoidance of contact with 
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wild birds. Upon vaccination, IBV generates antigenic vari¬ 
ants that can result in incomplete protection (Cavanagh 2005). 

Coronaviruses phylogenetically related to IBV were 
largely detected in several wild bird species in Northern 
England, as well as novel viruses that are not associated to 
classical IBV types (Hughes et al. 2009). This was also 
reported in the Beringian area between Alaska and Siberia 
(Muradrasoli et al. 2010) and in wild bird surveillance car¬ 
ried out in Hong Kong and Cambodia (Chu et al. 2011). 
These studies highlight the wide host range, geographic dis¬ 
tribution, and diversity of avian Co Vs, suggesting a possible 
impact on animal health and later a possible risk for human 
health (Muradrasoli et al. 2010). 

In the present work, we report the detection and the phylo¬ 
genetic aspects of coronaviruses detected in wild birds sam¬ 
pled at Alaotra Lake, Madagascar. This lake constitutes an 
exceptional site for the settlement of both endemic and migra¬ 
tory birds (BirdLife International 2014a). The sanitary and 
economic impacts of avian coronaviruses in the Madagascar 
poultry sector are still unknown. 

Material and methods 

For this work, we used part of the samples collected in the 
framework of the GRIPAVI project (http://gripavi.cirad.fr/ 
en/) launched in 2007, and concluded in 2011 by CIRAD 
in collaboration with the Departement de Recherche 
Zootechnique et Veterinaire du Centre National de la 
Recherche Appliquee au Developpement Rural (FOFIFA- 
DRZV) in Madagascar, with the support of the French 
Ministry of Foreign Affairs. Cloacal swabs collected in 
2011 from 357 free-ranging birds at Alaotra Lake on Mad¬ 
agascar highlands were used. The samples were dipped in a 
transport medium consisting of isotonic phosphate-buffered 
saline, pH 7.0-7.4, with antimicrobial additives (penicillin 
10,000 U/mL, streptomycin 10 mg/mL, amphotericin B 
25 pg/mL, and gentamycin 250 pg/mL) supplemented with 
20 % glycerol and stored in liquid nitrogen containers. 
They were shipped in dry ice and kept in the laboratory 
at -80 °C until processing. All samples were manipulated 
in a biosafety level 3 laboratory. 

Viral RNA was extracted from samples by a high- 
throughput-automated workstation Biomek FX P (Beckman 
Coulter) using the Nucleospin RNA virus kit, according to 
manufacturer’s instructions (Macherey Nagel). The viral 
RNA was resuspended in nuclease-free water and stored at 
-80 °C. Moloney murine leukemia virus reverse transcriptase 
(Invitrogen) and random hexamers were used in 20 pL reac¬ 
tions to generate cDNAs from 10 pL of the extracted RNA 
suspension, according to the manufacturer’s instructions. 

The obtained cDNAs were submitted to a pancoronavims 
nested PCR (nPCR) with Top Taq DNA polymerase (Qiagen), 


according to manufacturer’s instructions. A segment of the 
RNA-dependent RNA polymerase (RdRp) was targeted with 
the following primers: CoV-Fl (5'-GGKTGGGAYTAYC 
CKAARTG-3'), CoV-Rl (5 '-TGYTGTSWRCARAAYTCR 
TG-3'), CoV-F2 (5'-GGTTGGGACTATCCTAAGTGTGA- 
3'), and CoV-R2 (5'-CCATCATCAGATAGAATCATCATA- 
3') (Chu et al. 2011). Briefly, the nPCR was performed as 
follows: the first reaction was performed in a 25-pL volume 
containing 20 ng of cDNA with primers (CoV-Fl and CoV- 
Rl). The cycling conditions were 5 min at 95 °C, 40 cycles of 
45 s at 94 °C, 40 s at 49 °C, 45 s at 72 °C and a final incubation 
at 72 °C for 7 min. For the second (nested) reaction, we used 
2 pL of the first reaction under the same amplification condi¬ 
tions of the first PCR, but using primers CoV-F2 and CoV-R2 
in the assay. 

Standard precautions were taken to avoid PCR contamina¬ 
tion; blank controls without template were included in every 
set of five RT-PCR assays. PCR products were electropho- 
resed in 1.5 % agarose gels and the products visualized on 
UV light. The products with the expected size (400 to 
440 bp) were excised and purified using QiaQuick gel extrac¬ 
tion Kit (Qiagen) following manufacturer’s instructions. 

The positive samples were sequenced (Cogenics) with 
primers CoV-F2 and CoV-R2, using Big dye Terminator 
v3.1 (Applied Biosystems) cycle sequencing kit according to 
the manufacturer’s instructions. 

The BLAST was used to detect homologous regions in 
sequence databases. Sequences were aligned by ClustalW 
and MEGA version 6 (Tamura et al. 2013). Phylogenetic anal¬ 
ysis was performed based on a fragment containing approxi¬ 
mately 350 bp of the RdRp CoV gene using the neighbor¬ 
joining method with bootstrap support (1000 replicates) in 
MEGA version 6 (Tamura et al. 2013). 

Results and discussion 

Samples originating from 17 bird species were tested for 
Co Vs polymerase RdRp gene by RT-PCR, and 28 of the 
357 cloacal samples were positive (7.8 %). Positives were 
found in 11 species (Table 1). All sequences were submitted 
as GenBank:KM093872 to GenBank:KM093897; 
GenBank:KM093899 and GenBank:KM093901. Phylogenet¬ 
ic tree reconstructions enabled us to classify all the 
coronaviruses identified in this study as gammacoronaviruses 

(Fig. 1). 

These gammacoronaviruses were detected in a wide broad 
of bird orders, such as the following: Gruiformes, 
Passeriformes, Ciconiiformes, Anseriformes, and 
Charadriiformes. Nucleotide distances between Co Vs se¬ 
quences in this study were 0.0-23.5 %. Sequences formed five 
different subgroups within two main clusters, revealing the 
diversity among avian Co Vs strains from the same geographic 
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Table 1 Detection of gammacoronaviruses in avian samples from 
Alaotra Lake, Madagascar 


Order Species 

Sampled (n) 

Positive (n) 

Rate % 

Anseriformes 

Anas erythrorhyncha 

82 

4 

4.9 

Dendrocygna bicolor 

21 

0 

0 

Anas melleri 

15 

0 

0 

Dendrocygna viduata 

40 

3 

7.5 

Anas hottentota 

54 

1 

1.8 

Charadriiformes 

Gallinago macrodactyla 

18 

3 

16.6 

Charadrius pecuarius 

15 

3 

20 

Rostratula benghalensis 

4 

1 

25 

Ciconiiformes 

Bubulcus ibis 

8 

1 

12.5 

Ardeola ralloides 

3 

0 

0 

Cuculiformes 

Centropus toulou 

4 

0 

0 

Gruiformes 

Gallinula chloropus 

32 

3 

9.4 

Porphyrula alleni 

37 

5 

13.5 

Rallus madagascariensis 

8 

1 

12.5 

Dryolimnas cuvieri 

2 

0 

0 

Passeriformes 

Foudia madagascariensis 

11 

3 

27.2 

Strigiformes 

Tito alba 

3 

0 

0 

Total 7 17 

357 

28 

7.8 


area. Seven sequences from dabbling and whistling ducks 
formed a separate and distinct subcluster with a within group 
distance of 0.4-2.1 % and a mean distance between the sub¬ 
clusters of 9 %. Three sequences (KM093887 to KM093889) 
from waders and rail were found to be phylogenetically closer 
to various IBV strains. Notably, three sequences derived from 
weavers and one from small waterbird (KM093890 to 
KM093892 and KM093901) fell into a different cluster within 
other gammacoronaviruses related to geese in Russia 
(Muradrasoli et al. 2010) and a whistling duck in Cambodia 
(Chu et al. 201 1). These sequences are up to 24.6 % divergent 
on the nucleotide level in comparison to the sequences recov¬ 
ered and located in the other main cluster. 

In the present study, gammacoronaviruses were detected in 
7.8 % of the examined wild bird samples. This prevalence is 
higher by a magnitude of four than the one found in the same 
birds for avian Paramyxovirus type 1 (Cappelle et al. 2014). 
Sequences recovered from dabbling and whistling ducks were 
closely related to sequences isolated from gulls in Russia 
(Muradrasoli et al. 2010), sharing up to 93 % of similarity 
on the nucleotide level. The fact that some sequences were 
more closely related to IBV strains highlights the question on 


the original reservoirs of IBV, if domesticated or wild birds 
(Hughes et al. 2009). The majority of the sequences here de¬ 
tected were phylogenetically related to avian Co Vs detected in 
rock sandpipers ( Eroliaptilocnemis ). These related sequences 
were reported in the Commander Island, Russian Federation 
(Muradrasoli et al. 2010), highlighting not only its potential of 
geographic distribution but the interspecies transmission, as 
such, very similar viruses can be detected in different avian 
species: ducks, plovers, herons, weavers, waders, small water 
birds, and painted snipes. Noteworthy, this study shows that 
the same species can host genetically diverse 
gammacoronaviruses. 

It is remarkable the phylogenetic similarities among most 
of the sequences herein detected and others from Asia. The 
only possible explanation lies on the coastal migratory move¬ 
ment Africa-Madagascar. The direct route of migration of wa- 
terbirds to Madagascar is from the Wetlands in Eastern Africa 
(Dodman and Diagana 2006). Although it is hard to predict, it 
would be possible that many migratory species could have 
previous contact with migrating birds from Westem/Central 
Siberia, Balkans, Black Sea, and Central Asia which use the 
East Africa flyways to reach wintering areas to rest along the 
river systems crossing the Arabian Peninsula and the Nile 
(BirdLife International 2014b). 

Taken together, our results demonstrate that there is circu¬ 
lation of genetically divergent avian Co Vs among different 
wild bird species at Alaotra Lake, Madagascar. This lake is 
the island’s most important rice-growing region and a rich 
habitat for wildlife, including relatively many endemics, rare, 
and endangered species (Ferry et al. 2009). 

Due to the environmental conditions that prevailed in this 
region, it is assumed that in the case of introduction of other 
avian viruses, like avian influenza H5N1, this area may fa¬ 
cilitate virus spread and persistence (Guerrini et al. 2014). 
Considering our results on the detection of avian CoVs in 
the samples studied here, we might assume that these viruses 
are widespread among birds found in this environment. The 
detection of sequences related to gammacoronaviruses in 
this study area and in geographically distinct areas, such as 
Russia, Alaska, and Cambodia indeed indicates that Co Vs 
are widespread among birds associated with water environ¬ 
ments, raising questions on possible implications on wildlife 
and poultry production. 

The number of samples analyzed in this study was limited 
to predict the real prevalence of gammacoronavirus in Mada¬ 
gascar, but our results emphasize the circulation of genetically 
diverse coronaviruses in apparently healthy populations of 
wild birds from different species. Additional genetic data on 
coronavirus surveillance in a wide number of individuals and 
species, as well as the complete genomic characterization of 
these wild bird coronaviruses are required to better understand 
their evolution, the relationship between these viruses and 
wild bird population, and the risk of transmission to poultry. 
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Fig. 1 Phylogenetic relationship 
between the nucleotide sequences 
corresponding to a portion of 
350 bp of the Coronavirus gene 
(RNA-dependent RNA 
polymerase gene) amplified from 
wild birds at Alaotra Lake. 
Coronavirus detected in wild bird 
species in this study are denoted 
by black diamonds and their 
corresponding GenBank 
accession numbers. Sequences 
from others coronavimses were 
included for comparative 
purposes 
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